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Three dark matter direct detection experiments (DAMA/LIBRA, CoGeNT, and CRESST-II) have 
each reported signals which are not consistent with known backgrounds, but resemble that predicted 
for a dark matter particle with a mass of roughly ~10 GeV and an elastic scattering cross section 
with nucleons of ~10~''^-10~'"^ cm'^. In this article, we compare the signals of these experiments and 
discuss whether they can be explained by a single species of dark matter particle, without conflicting 
with the constraints of other experiments. We find that the spectrum of events reported by CoGeNT 
and CRESST-II are consistent with each other and with the constraints from CDMS-II, although 
some tension with xenon-based experiments remains. Similarly, the modulation signals reported 
by DAMA/LIBRA and CoGeNT appear to be compatible, although the corresponding amplitude 
of the observed modulations are a factor of at least a few higher than would be naively expected, 
based on the event spectra reported by CoGeNT and CRESST-II. This apparent discrepancy could 
potentially be resolved if tidal streams or other non-Maxwellian structures are present in the local 
distribution of dark matter. 
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I. INTRODUCTION 

Recently, a great deal of attention has been given 
to a number of dark matter experiments which have 
reported signals that do not appear to be consistent 
with any known backgrounds. The longest standing of 
these claims is from the DAMA (and more recently the 
DAMA/LIBRA) collaboration, which has reported an- 
nual variation in the event rate at their lowest observed 
energies, in keeping with the predictions from dark mat- 
ter scattering p]. Earlier this year, the CoGeNT collab- 
oration reported that their previously observed excess of 
low-energy events also exhibits seasonal variation [51 [3] , 
similar to the signal reported by DAMA/LIBRA. Most 
recently, the CRESST-II collaboration has reported an 
excess of events potentially attributable to dark mat- 
ter 0]. 

A casual comparison of these results can be confusing 
or even misleading. Conventionally, both experimental- 
ists and theorists in the field of dark matter direct detec- 
tion report their constraints and other results as derived 
using specific astrophysical assumptions and estimates 
of detector response. Adopting a single choice for each 
of these characteristics (even if that choice represents a 
reasonable estimate), rather than marginalizing over the 
possible or plausible range of those choices, can lead to 
regions of compatibility that are artificially small, and to 
constraints which are artificially stringent. 

In this paper, we revisit the signals reported by the 
CRESST-II, CoGeNT and DAMA/LIBRA collabora- 
tions, and attempt to determine whether they can be 
consistently explained with dark matter. In particular, 
in Sec. [ll] we study and directly compare the spectra 
of excess events reported by CoGeNT and CRESST-II, 
and from this conclude that (for reasonable astrophysi- 
cal assumptions) a dark matter particle with a mass of 



roughly 10-20 GeV and an elastic scattering cross section 
with nucleons of approximately (1 — 3) x 10"''^ cm^ could 
account for the excess events reported by each of these 
collaborations. We also note that a sizable fraction of 
this parameter space is consistent with the constraints 
placed by the CDMS-II collaboration. The constraints 
presented by the XENON-100 and XENON-10 collab- 
orations, however, remain in conflict, unless either the 
response of liquid xenon to very low-energy nuclear re- 
coils is lower than previously claimed, or the dark mat- 
ter's couplings to protons and neutrons destructively in- 
terfere for a xenon target. In Sec. |III[ we compare the an- 
nual modulation signals reported by the DAMA/LIBRA 
and CoGeNT collaborations. Again, we find good agree- 
ment between the results of these two experiments, but 
point out that under common assumptions (Maxwellian 
velocity distributions and velocity independent scatter- 
ing cross sections), the amplitude of the observed modu- 
lation requires the dark matter to possess a significantly 
larger (by a factor of approximately 3-10) elastic scatter- 
ing cross section than would be inferred from the spec- 
tra reported by CoGeNT and CRESST-II. In Sec. [iVAl 
we explore how this apparent discrepancy between the 
observed event rate and modulation amplitude could po- 
tentially be resolved by the presence of streams or other 
non-Maxwellian velocity structures in the local distribu- 
tion of dark matter, or by dark matter with a velocity- 
dependent scattering cross section with nuclei. In Sec.|Vj 
we summarize our results and draw conclusions. 



II. THE NUCLEAR RECOIL SPECTRUM 

The CoGeNT and CRESST-II collaborations have 
each reported an excess of nuclear recoil candidate events, 
difficult to attribute to known backgrounds. In this sec- 
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tion, we discuss the energy spectra of these events and 
the dark matter parameter space which could potentially 
account for these signals. 



A. CoGeNT's Event Spectrum 

CoGeNT is a P-type point contact germanium detec- 
tor with very low levels of electronic noise, enabling sen- 
sitivity to very low-energy nuclear recoils, and thus very 
low mass dark matter particles. Located in Northern 
Minnesota's Soudan Underground Laboratory, CoGeNT 
observes nuclear recoil events as ionization, and has thus 
far reported the results of 15 months of data collection, 
taken between December 2009 and March 2011 

In the upper-left frame of Fig. [l] we show the spec- 
trum of events reported by CoGeNT as a function of 
ionization energy (in keV-electron equivalent, keVee), af- 
ter subtracting the L-shell electron capture peaks (as de- 
scribed in Ref. [S]). Above about 1.5 keVee, the spectrum 
observed by CoGeNT is approximately flat and feature- 
less, and is thought to be dominated by Compton scatter- 
ing events. At lower energies, the observed rate climbs 
rapidly. While contamination from (non-rejected) sur- 
face events is expected to contribute significantly near 
threshold, it does not appear to be possible to account 
for the observed low-energy rate with this or other known 
backgrounds. 

To assess the hypothesis that the excess events re- 
ported by CoGeNT are the product of the elastic scatter- 
ing of dark matter particles, we compare CoGeNT's event 
spectrum to that predicted from dark matter. The spec- 
trum (in nuclear recoil energy) of dark matter induced 
elastic scattering events is given by [B] 
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where Nt is the number of target nuclei, moM is the 
mass of the dark matter particle, pdm is the local dark 
matter density (which we take to be 0.3 GeV/cm'^), v is 
the dark matter velocity in the frame of the Earth, 
is the velocity of the Earth with respect to the Galac- 
tic halo, and f{v,v'e) is the distribution function of dark 
matter particle velocities, which we will take to be the 
standard Maxwell-Boltzmann distribution for the time 
being: 
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The Earth's speed relative to the Galactic halo is given by 
Ve = UQ-t-WorbCOS 7 cos[a;(t— io)] where Vq = tio-l-12km/s, 
'^orb = 30km/s, cos 7 = 0.51, Iq is the date of the peak 
in the annual modulation (generally anticipated to lie 
within several weeks of late May or early June), and 
Lu — 27r/year. As a default choice, we adopt the com- 
monly used values of = 220 m/s and a Galactic escape 
velocity of 544km/s. This function should be thought of 



as a reasonable, but approximate, parametrization of the 
dark matter's true velocity distribution. Departures from 
a Maxwellian velocity distribution are in fact required for 
consistency with observed (and simulated) halo density 
profiles (see, for example, Ref. [5J and references therein). 
Such departures can non-negligibly impact the spectrum 
of dark matter induced events, and can significantly mod- 
ify the degree of seasonal variation in the rate [3 |H] • In 
Sec.lIV A| we will return to this issue and examine the ex- 



tent to which non-Maxwellian structures in the velocity 
distribution could potentially impact the signals reported 
by direct detection experiments. 

As the germanium isotopes which make up the Co- 
GeNT detector contain little net spin, we consider spin- 
independent interactions to generate the observed events. 
In this case, we have 
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where ^„ is the reduced mass of the dark matter parti- 
cle and nucleon (proton or neutron), an is the scattering 
cross section of the dark matter particle with neutrons, Z 
and A are the atomic and mass numbers of the nucleus, 
and fn.p are the couplings of the dark matter particle to 
neutrons and protons, respectively. Unless stated oth- 
erwise, our results are calculated under the assumption 
that fp = fn- The nuclear form factor, F{q), accounts for 
the finite momentum transfer in scattering events. In our 
calculations, we adopt the Helm form factor with param- 
eters as described in our previous work f5^. To convert 
from nuclear recoil energy to the measured ionization en- 
ergy, we have scaled the results by the quenching factor 
for germanium as described in Refs. [IQlllI] (Qce = 0.218 
at Eji =3 keV, and with the energy dependence predicted 
by the Lindhard theory 12J). 

Although the majority of surface events have been re- 
moved from the spectrum presented by the CoGeNT col- 
laboration through the application of a rise-time cut, this 
spectrum docs not take into account any inefficiencies in 
their surface event rejection algorithm. While it was ini- 
tially estimated that the signal would suffer only a mi- 
nor degree of contamination from non-rejected surface 
events [2], the CoGeNT collaboration has recently re- 
ported a somewhat higher estimate for the rate of non- 
rejected surface events near their energy threshold |13j . 
This (preliminary) estimate of the fraction of the event 
spectrum which consists of non-surface events {i.e. nu- 
clear recoil candidate events) is shown in Fig. [I Any 
unidentified surface events constitute an additional back- 
ground that should be accounted for when attempting to 
identify a dark matter signal from the CoGeNT detector. 
As can be observed from Fig. [2] this new measurement 
significantly reduces the estimate of the number of low- 
energy events that could potentially be attributed to dark 
matter. The precision of this measurement is expected 
to improve over time as more statistics are accumulated. 

In Fig. [T] we show the spectrum events at CoGeNT, 
for several different choices of the surface event correction 
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FIG. 1: In the upper-left frame, we show the raw spectrum of nuclear recoil candidate events as observed by CoGeNT, as 
originally presented in Ref. In the other three frames, this spectrum has been corrected using three different estimates for 
CoGeNT's surface event correcton factor, as shown in Fig. [2] In each frame, a spectrum of events from dark matter is shown 
(dashed line), along with this signal plus a flat background from Compton scattering (solid line). 
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FIG. 2: CoGeNT's surface event rejection correction factor (the fraction of nuclear recoil candidate events that are not surface 
events) as recently presented in Ref. [TH]. The four curves shown (including the horizontal line) correspond to the correction 
factors used to generate the corresponding spectra in Fig. [T] 



factor. As mentioned previously, the upper-left frame of 
Fig. [l] depicts the spectrum of events assuming a per- 
fect surface event rejection efficiency (100% of all surface 
events are identified as such, at all energies). The other 



three frames show the remaining spectrum of events after 
applying the most mild (green), central (red), and most 
stringent (blue) correction factor, as shown in Fig. [2j to 
the raw spectrum. For each of these choices of the surface 
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FIG. 3: The 90% (solid) and 99% (dashed) confidence level contours for the spectrum of events observed by CoGeNT, with 
each color matching the corresponding correction factors shown in Fig. [2] The CRESST contours (dot-dashed) denote the 95% 
confidence level regions. A dark matter particle with a mass of approximately 10-20 GeV and an elastic scattering cross section 
with nucleons of approximately (1 — 3) x 10"'*^ cm^ can account for the excess events reported by each of these experiments. 



event correction factor, wc find that the resulting spec- 
trum can be explained by an elastically scattering dark 
matter particle, although with slightly differing ranges 
of masses and cross sections. In each frame, we show 
an example of a good fit, with the dashed line denoting 
the signal from the dark matter alone and the solid line 
also including a flat background from Compton scatter- 
ing events, in each case correcting for the detector effi- 
ciency, as described in Ref. [2]. 

In Fig. [3] we plot the regions of dark matter parameter 
space that provide a good fit to the efficiency corrected 
spectra shown in Fig. [T] In each case, we have allowed 
the normalization of the fiat background to float. As 
expected, the inclusion of the non-rejected surface event 
background shifts the preferred region towards smaller 
cross sections, as there is now less dark matter signal then 
expected previously. The range of dark matter masses 
favored also shifts upward somewhat as a result of the 
additional surface event background. 

B. CRESST's Event Spectrum 

The CRESST-II collaboration makes use of eight 300 
gram cryogenic CaW04 detectors, operating in Italy's 
Laboratori Nazionali del Gran Sasso. Due to the rel- 
atively light oxygen and calcium nuclei in their target, 
CRESST is quite sensitive to dark matter particles in 
the mass range favored by CoGeNT. CRESST observes 
events through both scintillation and heat (phonons), 
enabling them to discriminate nuclear recoil candidate 
events from a variety of backgrounds. 

Very recently, the CRESST-II collaboration released 



an analysis of their first 730 kg-days of data, taken over a 
period between 2009 and 2011 p^. The analysis identified 
67 low-energy nuclear recoil candidate events, which is at 
least 30% more than can be accounted for with known 
backgrounds. The CRESST-II collaboration has assessed 
the statistical significance of this excess to be greater 
than 4(7. 

The CRESST-II analysis identified two distinct regions 
of dark matter parameter space which are compatible 
with the observed excess (see Fig. |3| . In the high mass 
region (referred to as Ml), the majority of the excess 
events arise from dark matter recoils with tungsten nu- 
clei. Within the low mass (M2) region, in contrast, the 
excess events are dominated by recoils on both oxygen 
and calcium. In an independent analysis based on the 
publicly available portions of the CRESST data, Ref. [13] 
identified a similar, but somewhat larger, region of com- 
patible dark matter parameter space. 

As can be seen in Fig. [3j the dark matter parameter 
space favored by CRESST-II is compatible with the re- 
gion implied by CoGeNT 's spectrum, after correcting for 
surface event contamination. In particular, a dark matter 
particle with a mass of roughly 10-20 GeV and an elastic 
scattering cross section with nucleons of (1 — 3) x 10"''^ 
cm^ could account for the excess events reported by both 
collaborations. 

C. Constraints From Other Experiments 

A number of direct detection experiments have pro- 
duced constraints which are relevant to the interpretation 
of the events reported by CoGeNT and CRESST-II. In 
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FIG. 4: A comparison of CoGeNT's spectrum (using the central estimate for the surface event correction factor, as shown in 
Fig. [2| with that reported by the GDMS-II collaboration [15]. We have subtracted the fiat (Compton scattering) component 
from CoGeNT's spectrum, and corrected for CoGeNT's efficiency. The solid curve represents the prediction for a 10 GeV dark 
matter particle with an elastic scattering cross section of a„ — 1.8 x 10^*^ cm^. The spectrum observed by the combination of 
all of CDMS's detectors is in good agreement with that observed by CoGeNT, although the spectrum from the single detector, 
T1Z5, is slightly lower than CoGeNT's below 1.2 keVee. 



particular, the impact of the constraints presented by the 
CDMS-II, XENON-100, and XENON-10 collaborations 
are significant for the regions of low-mass dark matter 
parameter space favored by CoGeNT and CRESST-II. 

The CDMS-II collaboration has presented the results 
of two analyses searching for dark matter particles in 
the mass range collectively favored by CoGeNT and 
CRESST [m [16]. Before taking into account the up- 
dated estimates of CoGeNT's surface event rejection ef- 
ficiency, these constraints appeared to be in conflict with 
a dark matter interpretation of CoGeNT's excess (see, 
however, Ref. 17\). As both CDMS and CoGeNT make 
use of germanium detectors, and thus are sensitive to sim- 
ilar systematic factors such as quenching factors for low- 
energy nuclear recoils, it was generally considered diffi- 
cult to reconcile CDMS's constraints with a dark matter 
interpretation of CoGeNT. In fight of the CoGeNT col- 
laboration recent estimate for their surface event rejec- 
tion efficiency, however, this apparent conflict seems to 
be largely resolved. In Fig. [4j we compare the spectrum 
at CoGeNT (after subtracting the flat, Compton scatter- 
ing, component, and applying the central estimate for the 
surface event correction factor) to that reported by the 
low-threshold analysis of CDMS-II. While the spectrum 
below 1.2 keVee from CDMS's T1Z5 detector is slightly 
lower than that observed by CoGeNT, the all-detectors 
spectrum reported by CDMS is in good agreement with 
CoGeNT's. 

The XENON-100 [Hj and XENON-10 [H] collabora- 
tions have also each reported rather strong constraints 
on the parameter space of low-mass dark matter parti- 
cles. As presented, these constraints appear to largely 
rule out the dark matter parameter space collectively 
favored by CoGeNT and CRESST. There are a num- 
ber of ways, however, in which these constraints could 
be significantly weaker than they might appear. Firstly, 



any uncertainties in the response of liquid xenon to very 
low-energy nuclear recoils (as encapsulated in the func- 
tions Leff and/or Qy) could significantly impact the cor- 
responding constraints for dark matter particles with a 
mass in the range of interest. The constraints from the 
XENON-100 collaboration were derived using measure- 
ments of the scintillation efficiency, Leff, as described in 
Refs. [20] , which have been criticized in Ref. [21] (see also 
Ref. [22)- Even modest changes to these values at the 
lowest measured energies (~3-4 keV) can lead to much 
weaker constraints on light dark matter particles. It has 
also been argued that the relatively large (9.3 eV) band- 
gap of xenon is expected to lead to a suppression of the 
response to nuclear recoils in the energy range of inter- 
est (see Ref. [32] and references therein). Many of these 
issues also apply to constraints on light dark matter mak- 
ing use of only the ionization signal in liquid xenon de- 
tectors [T^ . 

Alternatively, the constraints from XENON-100 and 
XENON-10 could be modified if dark matter particles 
do not have identical couplings to protons and neu- 
trons [33 |M| ■ In particular, for a ratio of couplings given 
by fnl fp ~ —0.7, the constraint from xenon-based ex- 
periments is weakened by a factor of ^-^20 relative to that 
found in the /„ = fp case ,24: . For this ratio of couplings, 
the cross section favored by CRESST-II would also be 
moved down by a factor of ~7 relative to that observed 
by CoGeNT. Alternatively, a ratio of fnl fp ~ —0.6 would 
reduce the strength of the XENON-100 and XENON-10 
constraints by a factor of 3-4, while also lowering the 
CRESST-II region (relative to that of CoGeNT) by a 
similar factor. 

Lastly, we note that a constraint has also been placed 
by making use of the CRESST commissioning run 
data [221 . These results appear to be in mild tension 
with the upper range (in cross section) of the parame- 
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ter space reported to be favored by the analysis of the 
CRESST-II cohaboration. 



III. ANNUAL MODULATION 

If a population of events observed in a detector are in 
fact the result of elastically scattering dark matter par- 
ticles, then we should expect the Earth's motion around 
the Sun to induce a degree of seasonal variation in the 
rate of those events. For most commonly assumed ve- 
locity distributions of dark matter particles, the rate is 
predicted to follow a roughly sinusoidal behavior, with 
a peak that occurs within several weeks of late May or 
early June pS] , 

The CoGeNT and DAM A/LIBRA collaborations have 
each reported the observation of annual modulation of 
their event rates. In this section, we characterize and 
compare the modulation signals reported by these col- 
laborations, and discuss whether these signals could be 
the result of dark matter. 



A. DAMA/LIBRA's Modulation 

The DAM A/LIBRA experiment makes use of a large 
mass detector (242.5 kg in its current form) consisting 
of high purity Nal(Tl) crystals, located at Gran Sasso. 
DAM A/LIBRA observes nuclear recoil events as scintil- 
lation, and is designed to search for time variations in 
their event rate, rather than to identify individual dark 
matter candidate events. 

Based on the data collected over a period of 13 an- 
nual cycles, the DAMA/LIBRA collaboration reports ev- 
idence of an annual modulation with a statistical signif- 
icance of 8.9(7. The variation of their rate is consistent 
with a sinusoid peaking at May 16±7 days at energies 
between 2 and 4 keV, May 22±7 days between 2 and 5 
keV, and May 26±7 days between 2 and 6 keV, consis- 
tent with that predicted for dark matter with a roughly 
Maxwellian velocity distribution. 

While DAMA/LIBRA's strategy of looking for an an- 
nual modulation in their rate can be successfully used to 
separate a dark matter signal from many possible back- 
grounds, one might worry about sources of background 
which could also exhibit seasonal variation. For exam- 
ple, the underground muon flux is known to modulate 
as a result of temperature variations in the stratosphere 
(although with a later phase and lower rate than is ob- 
served by DAMA/LIBRA (29j). Observed variations in 
the radon-induced background rate are also out-of-phase 
with the signal reported by DAMA/LIBRA. To date, no 
background has been identified with a phase, spectrum 
and rate compatible with DAMA/LIBRA's signal. 

The regions of dark matter parameter space in which 
the DAMA/LIBRA modulation can be accounted for 
depends strongly on the highly uncertain low-energy 
sodium quenching factor (the fraction of recoil energy 



of an elastic scattering event which is manifest as scintil- 
lation). In their analysis, the DAMA/LIBRA collabora- 
tion has adopted a canonical value of = 0.3 ± 0.01 
for this quantity, which they report to be the measured 
value averaged over the recoil energy range of 6.5 to 97 
keV [30]. Other groups have reported similar values: 
QNa = 0.25 ± 0.03 (over 20-80 keV), 0.275 ± 0.018 (over 
4-252 keV), and 0.4 ± 0.2 (over 5-100 keV) [H]. As the 
sodium quenching factor is generally anticipated to vary 
as a function of energy, it is very plausible that over the 
range of recoil energies relevant for light (5-20 GeV) dark 
matter particles (approximately 5 to 30 keV) the quench- 
ing factor could be quite different from the average values 
reported from these measurements. For recoil energies 
below approximately 20 keV, Ref. [35' reports a mea- 
surement of QNa = 0.33 ± 0.15, whereas Ref. [33j reports 
a somewhat smaller value of QNa — 0.252 ±0.064 near 10 
keV. The results of a very recent and preliminary mea- 
surement favor values of QNa ~ 0.15 — 0.2 at similarly low 
energies [131 [31] . At this time, we choose to keep an open 
mind regarding the relevant low-energy quenching factor 
for sodium, and will consider a range of values between 
QNa ~ 0.15 and 0.40. Based on theoretical considera- 
tions [35] |3S] and recently experimental evidence [T^ , we 
do not consider the possibility that channeling plays an 
important role at DAMA/LIBRA. 

For a quenching factor of QNa ~ 0.25 (0.15, 0.4), elas- 
tically scattering dark matter with a mass in the range 
of approximately 8-19 GeV (12-19 GeV, 6-14 GeV) can 
accommodate the spectrum of the modulation ampli- 
tude reported by DAMA/LIBRA (see Fig. [5|, assuming 
a Maxwellian velocity distribution with typical parame- 
ters (for earlier fits of DAMA data to light dark matter 
particles, see Refs. [371 133). The allowed regions do not 
extend to masses above about 18-20 GeV, where scatter- 
ing with iodine nuclei begins to dominate. Under these 
same assumptions, an elastic scattering cross section of 
CT„ ~ (0.7—3) X 10~^° cm^ is required to produce the ob- 
served magnitude of DAMA/LIBRA's modulation, which 
is significantly larger than the cross section implied by the 
spectra reported by CoGeNT and CRESST. Departures 
from a Maxwellian velocity distribution, however, could 
strongly impact (and potentially enhance) the observed 
modulation amplitude. We will return to this issue in 
Sec.|TV3 



B. CoGeNT's Modulation 

The CoGeNT collaboration has also reported evidence 
of an annual modulation in their event rate, although 
with a modest statistical significance of 2.8cr. Despite the 
lower statistical significance of this signal, it is interesting 
to compare the features of CoGeNT's time variation with 
that observed by DAMA/LIBRA. The peak of CoGeNT's 
phase is May 18±16 days, which is slightly earlier (at the 
1.6cr level) than that favored by DAMA/LIBRA. If the 
modulation signals reported by DAMA/LIBRA and Co- 
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FIG. 5: The 90% (solid) and 99% (dashed) confidence level contours for the spectrum of the amplitude of the annual modulation 
observed by DAMA/LIBRA, assuming a simple Maxwellian velocity distribution with vq = 220km/s and v^sc ~ 544km/s. 
Contours are shown for three choices of the low-energy sodium quenching factor, Qnh- 
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FIG. 6: A comparison between the modulation amplitude spectrum observe by DAMA/LIBRA and CoGeNT, independent of 
the dark matter's velocity distribution, following the approach of Ref. j27] . The comparison is done for a dark matter masses of 
10 and 15 GeV, and for three choices of the low-energy sodium quenching factor. The blue (red) error bars denote the CoGeNT 
modulation amplitude assuming a phase that peaks on April 18th (May 26th). The grey error bars denote the DAMA/LIBRA 
modulation spectrum. In normalizing the results, we have assumed the dark matter's elastic scattering cross section to scale 
with the square of the target's atomic number, A^. 



GeNT both arise from dark matter, a common phase that 
peaks in early May seems most likely [S] , in concordance 
with expectations for dark matter based on results from 
numerical simulations [7j. 

In comparing the spectra of the modulation amplitudes 



reported by DAMA/LIBRA and CoGeNT, it is possible 
to remove the dependence on the dark matter's velocity 
distribution, following the approach of Ref. (see also 
Ref. [2H|)- In Fig. [6| we present such a comparison. Al- 
though this comparison does not depend on the velocity 
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distribution of the dark matter particles, it does rely on 
assumptions pertaining to the mass of the dark matter 
particle, on the ratio of the elastic scattering cross sec- 
tions with germanium and sodium, and on the relevant 
quenching factors. Based on the shape of the CoGeNT 
and CRESST-II event spectra, we choose here to consider 
masses of 10 and 15 GeV, and assume a cross section 
which scales with of the target nucleus, as predicted 
for generic spin-independent scattering. In each of the 
three frames, we show the results for a different value of 
the low-energy sodium quenching factor. The spectrum 
of CoGeNT's modulation amplitude was determined us- 
ing the (publically available) 15 month CoGeNT data 
set [3], as described in Ref. |5]. 

From Fig. [6j it is immediately evident that the spec- 
trum and overall normalization of the modulation am- 
plitudes reported by CoGeNT and DAMA/LIBRA are 
quite similar. In fact, if the modulation reported by 
DAMA/LIBRA is the product of spin- independent elas- 
tic scattering with dark matter, then one should ex- 
pect CoGeNT to observe a modulation with broad fea- 
tures very much like that they report, and vice versa. 
The details of this comparison, however, depend signif- 
icantly on the value of the low-energy sodium quench- 
ing factor that is adopted. For a larger value (QNa ~ 
0.3 — 0.4) the steadily increasing modulation amplitude 
at low energies is seen by both experiments, while Co- 
GeNT's high modulation bin at ^2.5 keVee is not con- 
firmed by DAMA/LIBRA. As this feature is apparent 
only in one bin with a sizable error bar, we consider 
it possible that this bin represents a statistical fluctu- 
ation which may disappear with further data from Co- 
GeNT. Alternatively, if a lower sodium quenching factor 
is adopted (QNa ~ 0.15 — 0.2, as favored by Ref. p[3 l l34 ] ). 
the modulation reported by DAMA/LIBRA can overlap 
with CoGeNT's 2.5 keVee bin, while CoGeNT's lower 
energy modulation falls below the energy threshold of 
DAMA/LIBRA. In this case, one could consider the pos- 
sibility that this narrow feature results from a velocity of 
stream of dark matter present in the local halo. 



IV. WHY IS THE OBSERVED MODULATION 
AMPLITUDE SO LARGE? 

In the previous sections, we found that the event 
spectra observed by CoGeNT and CRESST-II are com- 
patible with arising from the same dark matter parti- 
cle. Similarly, the modulation amplitudes reported by 
DAMA/LIBRA and CoGeNT appear to be mutually con- 
sistent. Under the standard assumptions of a Maxwellian 
velocity distribution and velocity-independent scattering 
cross sections, however, the spectrum and rate of events 
reported by CoGeNT and CRESST-II would lead one to 
expect a signficantly smaller (by a factor of 3-10) modu- 
lation amplitude than is observed by DAMA/LIBRA and 
CoGeNT. In this section, we discuss how departures from 
these assumptions could explain why DAMA/LIBRA 



and CoGeNT have observed more modulation than would 
be naively predicted. 



A. Streams and Other Non-Maxwellian Velocity 
Distributions 

Numerical simulations of the formation and evolu- 
tion of Milky Way-like dark matter halos have become 
increasingly sophisticated in recent years. These sim- 
ulations find that although simple halo models with 
Maxwellian distributed velocities are likely to represent 
a reasonable zeroth order description of the distribution 
of dark matter in our galaxy, significant departures from 
such models are to be expected [55] . 

When considering relatively light dark matter parti- 
cles, as we are in this paper, the behavior of the velocity 
distribution near the escape velocity of the galaxy is of 
particular importance. In order to produce a measurable 
nuclear recoil, a low-mass dark matter particles requires 
greater speeds than would be necessary for a heavier par- 
ticle. As a result, a detector may only be sensitive to the 
high velocity tails of the dark matter distribution, where 
departures from Maxwellian behavior are expected to be 
most significant [5|. 

Small-scale structure of the Milky Way's halo can also 
play an important role in interpreting signals from di- 
rect detection experiments. The dark matter halo of our 
galaxy formed through a sequence of mergers of many 
smaller halos; a process known as hierarchial structure 
formation. High resolution simulations have found that 
many smaller halos survive this process and remain in- 
tact today, residing as substructures within larger ha- 
los [40]. Furthermore, many of these subhalos have a 
great deal of their outer mass stripped, resulting in the 
formation of cold tidal streams. The presence of such 
streams in or around the Solar System would introduce 
departures from the Maxwellian velocity distributions. 
While such streams could potentially effect the spectrum 
of dark matter-induced events that are observed in di- 
rect detection experiments [41], these effects are often 
far more pronounced in the modulation signals of such 
experiments [7]. The presence of such streams can sig- 
nificantly enhance a modulation signal, as well as shift 
the phase of the modulation relative to that predicted in 
more simple halo models [SI [T] . 

In Fig. [Tj we show examples of how streams might im- 
pact the spectrum of the modulation amplitude, as ob- 
served by CoGeNT. In the left frame, the results from 
a simple Maxwellian distribution are shown, for a dark 
matter mass of 10 GeV and an elastic scattering cross sec- 
tion of cr„ = 3 X 10"^-'^ cm^. As previously emphasized, 
the overall normalization of the predicted modulation is 
a factor of a few smaller than is observed. In the central 
frame, we add an additional tidal stream of dark mat- 
ter, with a velocity of 475 km/s, a local density of 0.072 
GeV/cm'^ (24% of the density of the smooth halo), and a 
dispersion of 15 km/s. The characteristics of this stream 
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FIG. 7: The impact of tidal streams on the modulation spectrum measured at CoGeNT. In the left frame, the result of a simple 
Maxwellian distribution is shown, for a dark matter mass of 10 GeV and an elastic scattering cross section of (t„ = 3 x 10~*^ 
cm^. In the center frame, we add a tidal stream of dark matter, with a velocity of 475 km/s, a local density of 0.072 GeV/cm^ 
(24% of the density of the smooth halo), and a dispersion of 15 km/s. The characteristics of this stream were chosen to explain 
CoGeNT's rather high error bar at 2.5 keVee (and, for the appropriate choice of the sodium quenching factor, the peak in the 
DAMA/LIBRA spectrum as well; see the upper right frame of Fig. 6|. In the right frame, we add a second stream, with a 
velocity of 165 km/s, a local density of 0.045 GeV/cm^ (15% of the density of the smooth halo), and a dispersion of 15 km/s. 



were chosen to explain CoGeNT's rather high error bar at 
2.5 keVee (and, for the appropriate choice of the sodium 
quenching factor, the peak in the DAMA/LIBRA spec- 
trum as well; see the upper right frame of Fig. [6]). In the 
right frame, we add a second stream, with a velocity of 
165 km/s, a local density of 0.045 GeV/cm^ (15% of the 
density of the smooth halo), and a dispersion of 15 km/s. 

Given the presently limited resolution of numerical 
simulations, it is difficult to assess the probability of sig- 
nificant tidal streams being present in our local neigh- 
borhood. Relatively small streams many order of mag- 
nitude below the length scales that can be currently re- 
solved could be very important. As an approximate lower 
limit, we note that current simulations _40J find signifi- 
cant streams to be present at our location of the Milky 
Way in roughly a few percent of realizations [7] . 



B. Velocity-Dependent Dark Matter Scattering 

If the dark matter's scattering cross section with nu- 
clei increases with the velocity of the dark matter parti- 
cle, the degree of seasonal variation in the observed rate 
can be larger than is predicted in the standard (velocity- 
independent) case. Inelastic dark matter scenarios are a 
well known example of models in which the dark mat- 
ter possesses velocity-dependent cross sections. In such 
models, the dark matter can only scatter with nuclei by 
being excited into a slightly heavier (typically on the or- 
der of 100 keV) state [32]. This requirement suppresses 
the rate of low energy events, and can increase the degree 
of annual modulation. 

Inelastic dark matter, however, does not appear help 
in reconciling the spectra observed by CoGeNT and 
CRESST-II with the modulation of DAMA/LIBRA and 
CoGeNT. In particular, the spectrum of events from in- 
elastically scattering dark matter is predicted to be quite 
flat at low energies, unlike that observed by CoGeNT and 
CRESST-II. 



Other models which introduce a velocity dependent 
scattering cross section include form factor dark mat- 
ter [43] and resonant dark matter models [44] (see also 
Refs. [45] [46]). Each of these classes of models hold 
promise for potentially explaining the large degree of 
modulation observed by DAMA/LIBRA and CoGeNT. 
In the first case, a new form factor is introduced which 
induces a momentum dependence in the interaction be- 
tween dark matter and nuclei, enhancing the cross sec- 
tion for higher velocity dark matter particles. While this 
feature can boost the observed modulation amplitude, it 
will also distort the spectrum of events. In the case of 
resonant dark matter, the interaction cross section is sig- 
nificantly enhanced near a particular center-of-mass en- 
ergy, leading to large (and potentially narrow in energy) 
modulation amplitudes. 

As the CoGeNT collaboration collects more data, 
the spectrum of the modulation amplitude will become 
rapidly better measured, making it possible to begin to 
discriminate between the various options described in 
this section. By the summer of 2012, CoGeNT will have 
doubled the size of its data set, and with less background 
contamination from L-shell electron capture peaks than 
was present in earlier data. In addition, the CoGeNT 
collaboration plans to deploy the first of four CoGeNT-4 
(C4) detectors in early 2012, roughly quadrupling their 
effective target mass. If streams or resonances are respon- 
sible for a significant fracton of the observed modulation, 
these features will become increasingly apparent as this 
data set grows. 

V. SUMMARY AND DISCUSSION 

In this article, we have compared the signals reported 
by the DAMA/LIBRA, CoGeNT, and CRESST-II col- 
laborations. We summarize our finding as follows: 

• The spectra of events reported by CoGeNT and 
CRESST-II are in good agreement, and (for a typi- 
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cal Maxwellian velocity distribution) are consistent 
with a dark matter particle with a mass of approx- 
imately 10-20 GeV and an elastic scattering cross 
section of cr„ « (1 — 3) x 10"'*"'^ cm^. This range 
of parameter space is roughly consistent with the 
constraints from CDMS-II, but is in tension with 
the constraints of xenon-based experiments unless 
the response of liquid xenon to very low-energy nu- 
clear recoils is lower than previously claimed, or the 
dark matter's couplings to protons and neutrons 
destructively interfere for a xenon target. 

• The spectra of the modulation amplitudes reported 
by DAM A/LIBRA and CoGeNT are also consistent 
with each other. Under the assumption of a typical 
Maxwellian velocity distribution, these modulation 
signals favor dark matter particles with masses of 
8-19 GeV and an elastic scattering cross section of 
a„ « (0.7-3) X 10-4° cm2. 

• The apparent mismatch between the elastic scat- 
tering cross sections required to produce the event 
spectra observed by CoGeNT and CRESST-II and 
those needed to produced the modulations reported 
by DAM A/LIBRA and CoGeNT could be poten- 
tially resolved if the local dark matter distribution 
contains streams or other highly non-Maxwcllian 
features, or if the dark matter's scattering cross 
section with nuclei is velocity-dependent. 

Taken together, these data appear to favor a dark mat- 
ter particle with a mass of approximately 10-15 GeV, 
and an elastic scattering cross section of roughly (t„ ^ 
2 X 10"'*^ cm^. This mass range is of particular interest 
in light of recent indirect detection results. In particu- 
lar, the spatial morphology and spectrum of gamma-rays 
observed from the Galactic Center can be explained by 
the annihilations of a 7-12 GeV dark matter particle, 
annihilating primarily to leptons, and with an annihila- 
tion cross section approximately equal to the value re- 
quired to generate the observed cosmological abundance 
of dark matter in the early Universe {av ~ 3 x 10"^^ 
cm'^/s) [47 , 48^. The same dark matter model (mass, an- 
nihilation cross section, annihilation channels, and halo 
profile) has also been shown to lead to the production of a 
diffuse haze of synchrotron emission consistent with that 
observed by WMAP It also appears that the ex- 

cess radio emission observed at higher galactic longitudes 
by the ARCADE 2 experiment [50] possesses a spectral 
shape and overall intensity consistent with originating 
from dark matter with the same mass, cross section, dom- 
inant channels, and distribution |511I52) . Lastly, we men- 
tion that '-^lO GeV dark matter particles with the same 
distribution and annihilation cross section would be capa- 
ble of depositing the required energetic electrons into the 
Milky Way's non-thermal radio filaments [53|, providing 
an explanation for their peculiar spectral features. Com- 
paring these resuhs to the CRESST-II, CoGeNT, and 
DAMA/LIBRA signals discussed in this paper, it may 



be the case that these experiments are each observing 
different facets of the same species of dark matter. 



Acknowledgements: While in the final stages of this 
project, Ref. [TJ appeared on the arXiv, which addresses 
many of the issues discussed in this paper. The authors 
are supported by the US Department of Energy. DH is 
also supported by NASA grant NAG5-10842. CK is sup- 
ported by a Fermilab Fellowship in Theoretical Physics. 



Bibliography 



[7 

[s: 

[lo; 
[11 

[12 

[13 
[14 
[15 

[16 

[17 
[18 

[19] 



R. Bernabei, P. Belli, F. Cappella et al, Eur. Phys. J. 
C6T, 39-49 (2010). arXiv:1002.1028 [astro-ph.GA]]. 
C. E. Aalseth et al. [CoGeNT Collaboration], Phys. 
Rev. Lett. 106, 131301 (2011). arXiv: 1002. 4703 [astro- 
ph.CO]]. 

C. E. Aalseth, P. S. Barbeau, J. Colaresi, J. I. Collar, 
J. Diaz Leon, J. E. Fast, N. Fields, T. W. Hossbach et al., 
Phys. Rev. Lett. 107, 141301 (2011). arXiv :1106.0650l 
[astro-ph.CO]]. 

G. Angloher, M. Bauer, I. Bavykina, A. Bento, C. Bucci, 

C. Ciemniak, G. Deuter, F. von Feilitzsch et al, 
[arXiv: 1109.0702] [astr o-ph.CO]]. 

D. Hooper, C. Kelso, [arXiv:11 06.l"066' [hep-ph]]. 

C. Savage, K. Freese, P. Gondolo, Phys. Rev. D74, 
043531 (2006). astro-ph/0607121 ; See also J. D. Lewin 
and P. F. Smith, Astropart. Phys. 6, 87-112 (1996). 
M. Kuhlen et al, JCAP 1002, 030 (2010) 
I arXiv: 09 12.2358 [astro-ph. G A] ] . 

M. Lisanti, L. E. Strigari, J. G. Wacker, R. H. Wechsler, 
Phys. Rev. D83, 023519 (2011). jarXiv: 1010.4300] [astro- 
ph.CO]]. 

C. Kelso, D. Hooper, JCAP 1102, 002 (2011). 
[arXiv: 101 1.3076 [hep-ph]]. 

For a summary of germanium quenching factor measure- 
ments, see S. T. Lin et al. [TEXONO Collaboration], 
Phys. Rev. D79, 061101 (2009). |arXiv:0712.1645 [hep- 
ex]]; Y. Messous, Astropart. Phys. 3, 361-366 (1995). 
P. S. Barbeau, J. 1. Collar and O. Tench, JCAP 09, 009 
(2007); P. S. Barbeau, J. I. Collar and P. M. Whaley, 
Nucl. Instr. Meth. A574, 385-391 (2007); P. S. Barbeau, 
PhD Thesis, University of Chicago (2009). 

D. Hooper, J. I. Collar, J. Hall, D. McKinsey, C. Kelso, 
Phys. Rev. D82, 123509 (2010). |arXiv: 1007. 1005 [hep- 
ph]]. 

Talk by J. Collar, TAUP 2011 Workshop, Munich, Ger- 
many, Sep. 5-9, 2011. 

J. Kopp, T. Schwetz, J. Zupan, [arXiv:1110.2721 [hep- 
ph]]. 

Z. Ahmed et al. [CDMS-Il Collaboration], Phys. Rev. 
Lett. 106, 131302 (2011). |arXiv:1011.2482 [astro- 
ph.CO]]. 

D. S. Akerib et al. [CDM S Collabora tion], Phys. Rev. 
D82, 122004 (2010). arXi v:1010.4290l [astro-ph.CO]]. 

J. I. Collar, arX iv: 1 1 03 . 34 8 1 ' [astro-ph. CO] ] . 

E. Aprile et al. [XENONIOO Collaboration], 

arXiv: 1104.2549 [astro-ph.CO]. 

J. Angle et al. [XENONIO Collaboration], Phys. Rev. 
Lett. 107, 051301 (2011). | |arXiv:1104.3688| [astro- 
ph.CO]]. 



BIBLIOGRAPHY 



11 



[20] G. Plante, E. Aprile, R. Budnik, B. Choi, K. -L. Gi- 
boni, L. W. Goetzke, R. F. Lang, K. E. Lim et al, 



arXiv: 1104.2587 [nud-ex]; 



J. Collar, arXiv: 1106.0653 [astro-ph.CO]. 

J. I. Collar, arXiv: 1010.5187 [astro-ph.IM]]. 

A. L. Fitzpatrick, D. Hooper, K. M. Zurek, Phys. Rev. 

D81, 115005 (2010). arXiv: 1003.0014 [hep-ph]]. 

J. L. Feng, J. Kumar, D. Marfatia, D. Sanford, Phys. 

Lett. B703, 124-127 (2011). arXiv:1102.4331| [hep-ph]]. 

A. Brown, S. Henry, IL Kraus, C. McCabe, 

'arXiv:1109.2589 [astro-ph.CO]]. 

A. K. Drukier, K. Freese, D. N. Spergel, Phys. Rev. D33, 
3495-3508 (1986). 

P. J. Fox, J. Liu and N. Weiner, Phys. Rev. D 83, 103514 
(2011) arXiv:1011.1915 [hep-ph]]. 

P. J. Fox, G. D . Kribs, T. M. P. Tait, Phys. Rev. D83, 

034007 (2011). larXiv:1011.19T0l [hep-ph]]. 

R. Bernabei, P. Belli, F. Cappella, R. CeruUi, C. J. Dai, 

A. d'Angelo, H. L. He, A. Incicchitti et al, AIP 

Conf. Proc. 1223, 50-59 (2010). arXiv:0912.0660 [astro- 

ph.GA]]. 

The DAMA CoUaboration, Nuovo Cim. 4, 26 (2003). 
See references contained within Ref. ^0, and K. Fushimi, 
H. Ejiri, H. Kinoshita et al, Phys. Rev. C47, 425-428 
(1993). 

D. R. Tovey, V. Kudryavtsev, M. Lehner et al, Phys. 
Lett. B433, 150-155 (1998). 

H. Chagani, P. Majewski, E. J. Daw et ai, JINST 3, 
P06003 (2008). arXiv:0806.1916 [physics.ins-det]]. 
J. Collar, private communication. 

N. Bozorgnia, G. B. Gelmini, P. Gondolo, JCAP 1011, 
019 (2010). arXiv: 1006.3110 [astro-ph.CO]]. 
C. Savage, G. Gelmini, P. Gondolo et al, Phys. Rev. 
D83, 055002 (2011). arXiv: 1006.0972 [astro-ph.CO]]. 
A. Bottino, F. Donato, N. Fornengo, S. Scopel, Phys. 
Rev. D69, 037302 (2004). |hep-ph/0307303 . 
P. Gondolo, G. Gelmini, Phys. Rev. D71, 123520 (2005). 
hep-ph/0504010. 

J. Diemand, B. Moore, J. Stadel, Mon. Not. Roy. As- 
tron. Soc. 352, 535 (2004). astro- ph/0402160. ; R. Wo- 
jtak, E. L. Lokas, S. Gottloeb^r^ "KT Mamon, 
_Mon. Not. Roy Astron. Soc. 361, LI (2005). f astro-" 
ph/0503391 ; S. H. Hansen, B. Moore, M. Zemp, 



J. Stadel, JCAP 0601, 014 (2006). astro- ph/0505420 ; 
S. H. Hansen, B. Moore, New Astron. 11, 333 (2006). 
|astro-ph/0411473] . 
[40] J. Diemand, M. Kuhlen, P. Madau, M. Zemp, B. Moore, 
D. Potter, J. Stadel, Nature 454, 735-738 (2008). 
(arXiv:0805.1244| [astro-ph]]; J. Diemand, M. Kuhlen, 



P. Madau, Astrophys. J. 657, 262-270 (2007). fastro- 
ph/061 1370; J. Stadel, D. Potter, B. Moore, J. Die- 
mand^ P. Madau, M. Zemp, M. Kuhlen, V. Quilis, 
[arXiv:0808.2981 [astro-ph]]; V. Springel, J. Wang, 
M. Vogelsberger, A. Ludlow, A. Jenkins, A. Helmi, 
J. F. Navarro, C. S. Frenk et al, Mon. Not. Roy. Astron. 
Soc. 391, 1685-1711 (2008). arXiv:0809.0898 [astro-ph]]; 
M. Zemp, J. Diemand, M. Kuhlen, P. Madau, B. Moore, 



[41] 



[42] 
[43] 

[44] 

[45] 

[46] 

[47] 
[48] 

[49] 



[50] 



[51] 

[52] 
[53] 



D. Potter, J. Stadel, L. Widrow, ,arXiv:0812.2033 [astro- 
ph]]. 

B. Moore, C. Calcaneo-Roldan, J. Stadel, T. R. Quinn, 
G. Lake, S. Ghigna, F. Governato, Phys. Rev. D64, 
063508 (2001). astro-ph/0106271 ; J. D. Vergados, 
S. H. Hansen, O Host, Phys. Rev. D77, 023509 (2008). 
(arXiv:0711.4895| [astro-ph]] . 

D. Tucker-Smith, N. Weiner, Phys. Rev. D64, 043502 
(2001). hep-ph/0101138; D. Tucker-Smith, N. Weiner, 
Phys. Rev. D72, 063509 (2005). |hep-ph/0402065 . 
B. Feldstein, A. L. Fitzpatrick, E. Katz, JCAP 1001, 020 
(2010) arXiv:0908.2991 [hep-ph]]; S. Chang, A. Pierce, 
N. Weiner, JCAP 1001, 006 (2010). arXiv:0908.3192| 
[hep-ph]]. 

Y. Bai, P. J. Fox, JHEP 0911, 052 (2009). 
|arXiv:0909.2900 [hep-ph]]. 

A. L. Fitzpatrick, K. M. Zurek, Phys. Rev. D82, 075004 

(2010) . arXiv: 10 07.5325 [hep-ph]]. 

M. Farina, D. Pappadopulo, A. Strumia, T. Volansky, 

I arXiv: 1107.0715 [hep-ph]]. 

D. Hooper, T. Linden, arXiv:1110.0006l [astro-ph.HE]]. 
D. Hooper, L. Goodenough, Phys. Lett. B697, 412-428 

(2011) arXiv:1010.2752 [hep-ph]]. 

D. Hooper, TT^Linden Phys. Rev. D83, 083517 
(2011), ,arXiv:1 011.4520' [astro-ph.HE]]; D. Hooper, 
D. P. Fi nkbeiner and G. Dobler, Phys. Rev. D 76, 083012 
(2007) arXiv:0705.3655 [astro-ph]]; D. P. Finkbeiner, 
arXiv:astro-ph/0409027 

D. J. Fixsen, A. Kogut, S. Levin, M. Limon, P. Lu- 
bin, P. Mirel, M. Seiffert, J. Singal et al, Astrophys. J. 
734, 5 (2011) arXiv: 0901.0555| [astro-ph.CO]]; M. Seif- 
fert, D. J. Fixsen, A. Kogut, S. M. Levin, M. Limon, 
P. M. Lubin, P. Mirel, J. Singal et al, Astrophys. J. 734, 
6 (2011) arXiv;0901.0559 [astro-ph.CO]]. 
N. Fornengo, R. Lineros, M. Regis, M. Taoso, 
iarXiv: 1108.0569 [hep-ph]]. 
M. Regis, private communication. 

T. Linden, D. Hooper, F. Yusef-Zadeh, Astrophys. J. 
741, 95 (2011) (arXiv: 1106. 5493| [astro-ph. HE] ] . 



